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We present the results of Monte Carlo simulations of two different Potts 
glass models with short range random interactions. In the first model a 
± J-distribution of the bonds is chosen, in the second model a Gaussian 
distribution. In both cases the first two moments of the distribution are 
chosen to be Jo — — 1) A J = +1, so that no ferromagnetic ordering of 
the Potts spins can occur. We find that for all temperatures investigated 
the spin glass susceptibility remains finite, the spin glass order parameter 
remains zero, and that the specific heat has only a smooth Schottky-like 
peak. These results can be understood quantitatively by considering small 
but independent clusters of spins. Hence we have evidence that there is 
no static phase transition at any nonzero temperature. Consistent with 
these findings, only very minor size effects are observed, which implies 
that all correlation lengths of the models remain very short. We also com- 
pute for both models the time auto-correlation function C (t) of the Potts 
spins. While in the Gaussian model C(t) shows a smooth uniform decay, 
the correlator for the ±J model has several distinct steps. These steps 
correspond to the breaking of bonds in small clusters of ferromagnetically 
coupled spins (dimers, trimers, etc.). The relaxation times follow sim- 
ple Arrhenius laws, with activation energies that are readily interpreted 
within the cluster picture, giving evidence that the system does not have 
a dynamic transition at a finite temperature. Hence we find that for 
the present models all the transitions known for the mean-field version of 
the model are completely wiped out. Finally we also determine the time 
auto-correlation functions of individual spins, and show that the system 
is dynamically very heterogeneous. 

PACS numbers: 64.70.Pf, 75.10.Nr, 75.50.Lk, 02.70.Lq 
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INTRODUCTION 



Recently much attention has been given to the study of several models of mean 
field spin glasses [p], §, |3|, |), such as the spin glasses with p— spin interaction || § 
and p-state Potts spin glasses with p > 4 0, 0, §, g, [T^, 0, [13], 0, P, ||, 



17| , [L8| , plgf] . The interest in these latter models is related to the fact that they 
exhibit two distinct transitions: A dynamic transition from ergodic to non-ergodic 
behavior at a temperature To where the phase space of the systems splits into 
many valleys separated by infinitely high barriers of the free energy , and a static 
transition at a temperature T < Tjj. Although at the static transition a (static) 
glass order parameter appears discontinuously and the internal energy and entropy 
have kink singularities, there is no latent heat (Fig. [I]). Due to these properties it 
has been proposed that these spin glass models are generic models for the structural 
glass transition from a supercooled fluid to an amorphous solid [^Cj, |2l], |23|, |24], 
25]. In fact, the slowing down of the relaxation near To is described by the same 
type of equations f|, 0, 0, |15|] as has been proposed before in the framework of 
the so-called idealized mode coupling theory (MCT) |26|, |27|, |28], to describe 



the onset of the slowing down in supercooled fluids near the glass transition, and 
it is known that MCT describes a variety of experiments and simulations quite 
successfully 
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In addition, the T— dependence of the entropy of these 
models (Fig. [I]) is reminiscent of the behavior of the entropy of supercooled fluids 
near the glass transition, which can be extrapolated to the Kauzmann temperature 
Tjc (which is below the experimental glass transition temperature T g ) []3T| at which 



it seems to be equal to the entropy of the corresponding crystal. Since one does 
not expect that the entropy of the (disordered) fluid is less than the entropy of the 
(ordered) crystal, one often invokes a static phase transition from the (metastable) 
supercooled fluid to a (metastable) ideal glass at a temperature T > T K p2l. |33|. 



However, the existing theories for such a static transition underlying the 
structural glass transition rely on mean-field type approximations whose validity 
is questionable [BS, 36]. Since real systems have short range interactions, it is 
relevant to ask which features of the mean field Potts glass (the scenario of Fig. [I] 
holds true for infinite range interactions) can also be found in systems with short 
range interactions. In fact, one must expect that in a short range system the 
temperature Tjj can only have the character of a crossover temperature, where 
valleys in phase space appear that are separated by barriers of large but finite 
height. Such a scenario is also suggested by the extended form of MCT |28f and 



is certainly also required from an analysis of experimental data |29[ which clearly 
rules out a true power-law divergence of the relaxation time at T D in supercooled 
fluids. Instead one finds in experiments a smooth crossover from an apparent 
power-law for T > To to an activated behavior (Arrhenius-behavior: lnr cx 1/T) 
for T < Td or an even stronger T— dependence. Of course, if the singularities 
at Tp are rounded off in systems with short range interactions it is of interest to 
ask whether the static transition at T is still observed or is also wiped out by 
the fluctuations that are neglected in mean field theories. Eastwood and Wolynes 



have recently put forward some arguments by which supposedly it can be decided 
whether or not a short range model shows close to T the behavior expected from 
the mean field theories |J7). However, these arguments rely on the idea that it is 
possible to partition a given configuration into droplets which have a well-defined 
free energy and surface tension. To what extend such ideas are indeed applicable 
also to simple spin glasses with nearest neighbor interactions is presently an open 
question, however. 

The motivation for the present study is to address some of the issues discussed 
in the previous paragraph by means of numerical simulations. The systems con- 
sidered are d = 3 dimensional 10— state Potts glasses with nearest neighbor inter- 
actions p8 |. In Sec. |T| we will describe the models that are studied in the present 



paper and define the quantities that are calculated. The main effort of our work 
is concentrated on the ± J model, and the results on the static and dynamic prop- 
erties of that model will be presented in Sees. |T| and fV[ Finally, in Sec. |V| we 
will present some results on the Gaussian model, while in Sec. |VI] we summarize 
our conclusions. 



II. MODELS AND DETAILS ON THE SIMULATION 

In this section we define the systems that are studied in the present work 
and introduce the observables that we will consider. Also some details on our 
simulation procedures are given. 

The Hamiltonian of the p-state Potts glass is given by 

w = -E JiAv^-i) , (i) 

where the "exchange constants" Jy are quenched random variables, and the symbol 
means that each pair of nearest neighbor sites i, j on the simple cubic lattice 
is included in the sum only once. In the present work we will always use p — 10, 
and will set the Botltzmann constant ks — 1. Most of our work concentrates on a 
bimodal distribution 

P(Jij) = x5(Jij - J) + (1 - x)5(Jij + J) , (2) 

since in a simulation of this model the transition probabilities used in the Monte- 
Carlo algorithm for the assignment of a new value of a Potts spin (we used a single 
spin-flip type algorithm) can be pre-computed and stored in a small table. Hence 
they need not be recomputed in the course of the simulation, which in turn leads 
to a faster algorithm. 

If we denote the first two moments of the distribution P(Jij) as J an d (AJ) 2 , 
respectively, 

J = [J«U (AJ) 2 = [4] av -([J tJ ] av ) 2 , (3) 



we obtain the relation with x and J from Eq. (|2]) as 
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x = -(l + J /J) , J=^/J$ + (Ajy . (4) 

(Note that the brackets [.] av in Eq. (||]) stand for the average over the distribution 
P(Jij) from Eq. (§). In practice this was done by using typically 100 independent 
realizations of bond configurations.) 

To fix the temperature scale we set A J = 1 and furthermore we choose Jo = — 1 
to avoid the tendency to ferromagnetic ordering. Note that in the Ising spin 
glass [0] the choice J = is most natural and does not lead to any ferromagnetic 
correlations in that model. However, this is no longer true in Potts glasses with 
p > 3 as we will show now: Consider the ferromagnetic susceptibility 

Xfm = iV[(m 2 )] av , (5) 

where m is the "magnetization" per spin of the Potts model (using the simplex 
representation discussed below) and iV = L 3 is the number of sites in the Lx Lx L 
simple cubic lattice. (Note that in order to minimize finite size effects we always 
used periodic boundary conditions.) The brackets (.) stand for the usual canonical 
average which were calculated by standard Monte Carlo averaging |J9[ |40| , using 
runs of a typical length of 10 8 Monte Carlo steps (MCS) per site. (Before these 
production runs were started, we equilibrated the samples for the equal amount 
of time.) 

The temperature dependence of xfm for the case J = and J = —1 is 
shown in Figs. ^|a and 0b, respectively. For J = 0, xfm gets rather large at 
low temperatures and increases rapidly with increasing L. In contrast to this the 
curves for Jo = — 1 stay on the order of unity in the whole temperature range and 
show, within the statistical errors, no systematic dependence on the lattice size. 
Hence we conclude that the choice Jo = — 1 does not lead to any ferromagnetic 
ordering of the system. 

According to Eq. (f|) the choice A J = 1 and J = — 1 leads to 

x=(l- l/V2)/2 w 0.146 and J = y/2, (6) 

and these are the values that we use in this work. 

In addition to the ±J model defined by Eqs. (|1|) and (0), we also investigate 
the Gaussian model defined by 



p{Jli> = Tsko exp 



(Jij - JoY 



(7) 



2(AJ) 2 

Results for this model will be discussed in Sec. [V|. 

The static quantities determined in our simulations include the energy per spin 
and the specific heat calculated from energy fluctuations, 

e(T) = ±[(H)U and c(T) = { [(tt 2 )] av - [(H}\ v } . (8) 



For defining observables like the magnetization, the glass order parameter, and 
the time-dependent spin auto-correlation function, etc., it is useful to choose a 
representation for the spins that takes into account the symmetry between their 
p possible states. This can be achieved by the so-called "simplex representation" 
41], |42]1 , in which each of the p states corresponds to a (p — 1) -dimensional vector 



S\ pointing towards the corner of a p— simplex, i.e. 

S x -Sy = (p5 xx/ -1) with A, A' e{l,2,...,p}. (9) 

To define a spin glass order parameter, we follow the standard method used for 
Potts glasses [17], |18|, |l9[ |43| by considering two real replicas a and (3 of the 
model system, i.e. two systems that have identical bond configurations, and to 
carry out for each of them an independent Monte Carlo simulation. The order 
parameter tensor is then defined as 



1 - 

-j2(s ha ns hP r we {1,2, 
i=i 



,p-i} 



(10) 



In an equilibrium simulation of a finite system with no external fields that couple 
to components of the order parameter, there is no symmetry breaking, and hence 



it is sufficient to consider the root mean square order parameter [IB, [IT], [18|, |T^, (43f : 



P 



1 p_1 
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By calculating a histogram of q (i.e. one always takes first the thermal average 
and then the sample average over the disorder) we can estimate the distribution 
P(q) of the order parameter. The various moments of this distribution are then 
given by 



i(q k )U = J q k P(q)dq . 



(12) 



The spin glass susceptibility can now be defined in terms of the second moment 

N 



Xsg 



p — 1 



[{q% 



(13) 



the reduced fourth-order cumulant |L6l 1/, 18, 19, 



Further quantities that exhibit interesting behavior if a phase transition occurs are 

(14) 



[(<!% 



(p-l)* [<g 2 >] 



u 



and the so-called Guerra parameter 



gwt)= [ S:-IS • (i5) 



J av 



Both quantities are defined such that for N —> oo they are zero throughout the 
disordered phase, while they differ from zero in the ordered phase. If the system 
has a phase transition at a nonzero temperature T c , one expects that the curves 
g 4 (iV, T), or G(N, T), plotted vs. T for different choices of N intersect at a common 
intersection point at T = T c . In particular, G(N,T) is a measure for the lack of 
self- averaging p| . 

To study the dynamical properties of the system, we calculate the time auto- 
correlation function of the Potts spins, 

1 N 

c{t) = N(p-i) D^o-^+^u • (16) 

(Note that in equilibrium, the right hand side of this equation depends only on 
the time difference t.) 

It is also of interest to study a time correlation function that is analogous to 
C(t), in which, however, one considers only the time auto-correlation function of 
the spin for each lattice site i, since such a quantity will yield information on the 



dynamical heterogeneity jy| [47], g§] of our model system. Below we will investigate 
these functions C»(t) in detail. 

To propagate the system in configuration space, we used the heat bath algorithm 



3£j. I. e., to every possible state i of a Potts spin (£ e {1, 2, • • • ,p}) a probability 
Pi is assigned, 

exp(-E e /T) 

Pi = — (17) 

E exp(-E m /T) 

m=l 

where Ei is the total energy of the system if the considered Potts spin is in state 
i. At each time step a random spin i is picked, the probabilities given by Eq. ( |i~7| ) 
are determined, and the spin is flipped to state £ with probability pg. (Hence the 
algorithm is quite efficient in finding energetically desirable states even for large 
values of p, in contrast to the standard Metropolis algorithm.) Note that in the 
present case the computation of Ei involves only a consideration of the z = 6 
nearest neighbors, unlike the mean field model studied previously 1TJ|, O, |18L ff9|. 



III. STATIC PROPERTIES OF THE ±J POTTS GLASS 

In Figs. | and £| we show the temperature dependence of the energy per spin 
e(T) and the specific heat per spin c(T), Eq. (|8|). We see that for T < 2 the energy 
seems to become independent of T, i.e. it must be close to the ground state energy 



e(0) ~ —8.8. This value for e(0) is not far from the value that one would predict 
if only the energetically favored ferromagnetic bonds were fully satisfied. In fact, 
in such a case we would have (note the normalization of "H in Eq. (|10D) 



e(0) = -3[Jx(p - 1) - J(l - x)(-l)} = -3J[xp - 2x + 1] « -9.2. (18) 

This estimate is only 5% lower than the actually observed value, showing that 
the effect of frustration in the present model is very small, unlike the case of the 
±J Ising spin glass [p]]. Since this frustation is so low, it should be possible to 
describe the thermodynamic behavior of the system as a superposition of small 
independent clusters of spins. In the Appendix we present such an approach and 
give the T— dependence of the energy and the specific heat (Eqs. (|A.3|) , (|A.4|) , and 
( |A.6| )). The result of this calculation is included in Figs. |3| and £| as well. We see 



that for the energy the theoretical curve at low T is somewhat lower than the data 
from the simulations (since not all effects of the frustration are taken into account) 
but that for intermediate and high temperatures the agreement is very good. 

Figure f| shows that for the case of the the specific heat the agreement between 
the data of the simulation and the theoretical calculation is even better. Note that 
the shape of c(T) resembles the one expected for a Schottky anomaly and a closer 
look at Eq. ( |A.4j ) shows that the full theoretical expression is indeed very similar 
to such a functional form (which is c(T) oc (3 2 exp(— pJ/T)). It is found that the 
dominant contribution to c(T) does come from clusters in which exactly two spins 
are coupled ferromagnetically and hence for which the energy needed to flip one 
spin is pJ. 

Note that using the theoretical approach presented in the Appendix it is also 
possible to calculate the ferromagnetic susceptibility shown in Fig. [|b. We have 
found, however, that the theoretical predictions are only in qualitative agreement 
with the results from the simulations and that significant quantitative differences 
exist. The reason for this quantitative discrepancy seems to be that this type of 
approach is not reliable in the presence of anti-ferromagnetic couplings, as a similar 
calculation for a pure anti-ferromagnetic model (with no disorder) has shown. 

Note that from Fig. |] and [| we can also conclude that within the accuracy of 
the data the energy as well as the specific heat are independent of the system size, 



in contrast to our findings for the mean-field version of the same model [18|, [19 
Hence we can conclude tentatively that the system does not show a phase transition 
at a finite temperature. 

The same conclusion emerges from the moments of the spin glass order param- 
eter distribution (Figs. ^|and []). In particular we see that the first moment shows 
a trivial iV -1 / 2 dependence and that xsg shows no sign of a divergence in the 
whole temperature regime investigated. At the lowest temperatures xsg seems to 
become constant, which is evidence that the relevant T— range has been probed. 
The value of the constant for T — > is around 3.0. Since for a short range system 
Xsg(^) can be written as a sum of the spin glass pair correlation function over 
all distances, this result again implies that these correlations are small and very 
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short ranged. For high temperatures the T— dependence of xsg is compatible with 
a 1/T— dependence (see inset of Fig. ^j, a functional form that is expected for a 
system below its critical dimension j49|. Also remarkable is the complete absence 
of finite size effects in this quantity, in strong contrast to the mean field version of 
the model 0,0. 

Finally we show in Figs. [?| and |] the temperature dependence of the cumulant 
g A {N,T) {Eq. QIJ)} and the Guerra parameter G(N,T) {Eq. (^)} for different 
system sizes. One recognizes that with increasing L these quantities converge to 
zero for all temperatures, which is again consistent with the complete absence of 
a finite transition temperature for this model. 



IV. DYNAMIC PROPERTIES OF THE ±J POTTS GLASS 



Having discussed the static properties of the ± J Potts glass, we now investigate 
its dynamical properties. 

In Fig. [5] we show the time dependence of the (average) spin-autocorrelation 
function C(t), see Eq. (|16|), for all temperatures investigated. Although the data 
shown are for L = 10, we found that within the accuracy of our data the curves 
are identical to the ones for larger system sizes. This can be seen in the inset of 
the figure where we compare for T = 1.6 the curve for L — 10 with the one for 
L = 16. Note that this absence of finite size effects is in strong contrast to the 
results found for the mean field version of the present model |1| 19 1. 



From the figure we see that the time-dependence of C(t) is rather unusual in 
that at low T the function does not show a two step relaxation process as in Ising 
spin glasses, but instead shows two plateaus: One with a height of approximately 
0.6 and another one with height 0.1. Note that normally the presence of such a 
plateau is associated with the existence of some sort of temporary cage in which 
the spins are trapped for a certain time. Thus the fact that we now see two 
plateaus might be taken as an indication that there are (at least) two different 
type of cages. However, as we will show below, for the present model the reason 
for the existence of the plateaus is rather different, despite their similarity with the 
plateaus seen in the mean field version of the same model close to the dynamical 



critical temperature Tjj [pi [TO, O, ELR |I9||. Last not least we remark that the 



existence of further plateaus cannot be ruled out since within the accuracy of our 
statistics and within the temperature range accessible in our simulation it is not 
possible to make any reliable statements on this. 

In order to measure the lifetimes of these plateaus, we define relaxation times 
Ti via the conditions 



C(t = n) = 0.4, C(t = r 2 ) = 0.08, and C(t = r 3 ) = 0.05. (19) 

Thus, the relaxation times Tj are obtained from the intersections of the relaxation 
functions C (t) at the various temperatures with the horizontal dashed straight lines 
shown in Fig. ^. The time T\ measures the lifetime of the first plateau, whereas r 2 
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and T3 measure the one of the second plateau. (The reason for using two different 
definitions for the latter will be discussed below.) 

Many glass forming systems have the property that their time correlation func- 
tions obey the so-called time-temperature superposition principle (TTSP) [j27f . 
This means that if one plots a time correlation function versus t/r, where r is 
the a— relaxation time, one finds that the correlators for different temperatures 
collapse onto a master curve. The existence of the TTSP is one of the main pre- 
dictions of the mode-coupling theory of the glass transition and it follows directly 
from the solution of the mode-coupling equations |27J . Since the latter are identi- 
cal to the one derived for the time-dependence for C(t) for the mean field version 
of the present model, it is of course of interest to check to what extend the TTSP 
holds for the present short range model as well. That this is indeed the case is 
demonstrated in Fig. [H] where we plot C{t) vs. t/ri and vs. t/r 2 . We see that 
with this type of plot the curves for the different temperatures do indeed fall on 
top of each other, as expected from the TTSP. (We mention that a plot of C{t) 
vs. t/r 3 also leads to a superposition of the curves.) 



In Fig. |TT] we show an Arrhenius plot of the temperature dependence of the 
relaxation times Tj. From this plot we see that the T— dependence of the times r 2 
and T3 is the same and hence we conclude that the exact definition how we defined 
these relaxation times is irrelevant. 

Also included in the graph are straight lines which represent fits to the 
data. Since these fits give a very good description of Tj we conclude that the 
T— dependence of Xj is just an Arrhenius law: 

Ti oc exp{E$/T) . (20) 

The activation energies are m 14.6 and E^ = E^' m 28.2. The interpretation 
of this observation is in fact quite simple: From Eqs. ([[]) and (§) we conclude that 
the breaking of one ferromagnetic bond costs an energy py/2 « 14.14 « E { 2\ 
while overturning a Potts spin which is coupled by two ferromagnetic bonds to 
neighboring spins would require an energy 2p^2 « 28.3 and this is Ef = Ef. 

In order to understand this behavior in more detail, it is illuminating to ana- 
lyze the relaxation dynamics of single spins. Note that for such an analysis it is 
necessary to simulate the system for a time that is on the order of 10 3 larger than 
the longest relaxation time of the system, since only then it is possible to obtain 
a statistics that is sufficiently accurate. Due to this heavy numerical task, we will 
present the following analysis only for one realization of the disorder. However, 
the results for two different samples gave qualitatively the same results. 



In Fig. [12] we show the time dependence of Cj(i) for three different temperatures 
(T = 2.4,2.0 and 1.6). Since L = 10, each plot has 10 3 different curves. One sees 
immediately that the relaxation dynamics of the system is very heterogeneous in 
that the typical decay time of the different curves spans many decades and that 
this range increases rapidly with decreasing T. We find that about 40% of the 
curves relax to (basically) zero very quickly, i.e. within just a few MCS. This fast 
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relaxation is the reason for the presence of the first plateau seen in C(t) at a height 
of approximately 0.6 (= 1 — 0.4). For T = 2.4 (Fig. |T2|a) most other curves decay 
to zero within a time of 10 3 MCS. However, if the temperature is decreased to 
T = 1.6, see Fig. |l|c, we see that the typical relaxation time of these curves has 
increased to about 10 4 MCS. Thus we conclude that the second step in C(t) is due 
to the relaxation of these spins. In addition we also recognize from the figure that 
some spins relax only on the time scale of 10 6 — 10 7 MCS. Thus the existence of the 
second plateau in C (t) is related to the presence of these slow spins. We also note 
that the presence of dynamical heterogeneities has also already been observed in a 
short range Ising spin glass JsTJ, [5l|. However, in that case no simple interpretation 
of the results has been given so far. 

Apart from the curves Ci(t) that, at low T, can be grouped together in a simple 
way, there are also curves which show a rather complicated time dependence in that 
they may show a fast decay at short times, one or two plateau(s) at intermediate 
times, before they decay to zero at long times. Since the total fraction of such 
curves is only a few (3-4) percent, and the number of such curves that relax only 
on the longest time scales does not seem to increase with decreasing T, we feel 
that they are not important for the average relaxation dynamics of the system. 
Therefore we will neglect them in the following. 

We have found that the shape of the remaining curves can be described very 
well by a simple exponential. In addition the T— dependence of relaxation time r[ 
of each curve is close to an Arrhenius law and therefore we can make the following 
Ansatz: 



t[ = a t expfo/T), ie{l,...,JV} • (21) 

A histogram with the distribution of the parameters a« and hi is shown in Fig. |13"| . 
Note that since we have determined the r[ at three different temperatures we can 
use the combinations (T = 1.6 and T = 2.0) as well as (T = 2.0 and T = 2.4) to 
determine the two parameters a« and fej. However, as can be seen from Fig. [13], 
the two combinations of temperatures give rise to very similar histograms which 
is evidence that the Ansatz fl2~l~|) holds and that the way the parameters a« and 6« 
were determined is irrelevant. 

While the distribution of a« seems to be continuous without any special features 
(see inset), the one for hi shows clearly three peaks: One at 6j ~ 0, one at at pJ 
and a third one at 2pJ respectively. Thus we see that the three different types of 
groups of spins identified in Fig. can also be found in this analysis. 

It turns out that all these features of the relaxation behavior can even be un- 
derstood quantitatively. First we mention our observation that all of the spins that 
relax quickly do not have ferromagnetic bonds. Recalling that each spin has z = 6 
neighbors and that there are p = 10 states at disposal, the probability that a spin 
has only antiferromagnetically coupled neighbors is given by (1 — s) 6 « 0.39, which 
accounts for the numerically estimated fraction of 40% of fast relaxing spins. More 
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general, the probability to have exactly k ferromagnetic bonds out of n bonds is 

p ^=w^w^ {1 - x) ■ (22) 

So, a first estimate would be to say that a spin with exactly k ferromagnetic bonds 
relaxes with a barrier given by kpj. Noting that P(l,6) ~ 0.40, P(2,6) ~ 0.17, 
P(3, 6) ~ 0.04 etc., we note that this approach would predict a second plateau 
at a height of 0.21 (resulting from the sum of probabilities with k > 2), and also 
predict a third plateau with a rather large height of about 0.046. Since we do not 
observe such a third plateau it is necessary to refine the approach: While spins 
with a single ferromagnetic bond need to overcome a barrier of pJ, the spins with 
two ferromagnetic bonds do not necessarily need to overcome a barrier of 2pJ. 
E.g., if one of the two spins that are connected ferromagnetically has itself only 
one ferromagnetic bond, it will relax with a barrier pJ, and thus can subsequently 
permit to the first spin with two ferromagnetic bonds also to relax, overcoming a 
barrier of p J only. The probability to have such an arrangement of bonds is given 
by P(2,6), that is the probability for the spin to have exactly 2 ferromagnetic 
bonds with its neighbors, times the probability that one of these two spins has 
no further ferromagnetic bond with its remaining 5 neighbors, and which is thus 
P(0, 5). This yields P(2, 6) • P(0, 5) ~ 0.077, so that the total percentage of spins 
relaxing with barriers higher than pJ is lowered to 0.13, very close to the value 
found for the second plateau from Fig. ||. Extending this reasoning also to spins 
relaxing with a barrier 3pJ, one estimates that the height of a third plateau would 
be around 0.006, and this is clearly too small to be observable numerically within 
the accuracy of our data. Note also that these numerical estimates completely 
ignore the small effects due to the frustrated bonds, which still must be present, 
as the estimate for the ground state energy shows. Nevertheless, these arguments 
show that the relaxation of the model is to a large extend accounted for by the 
behavior of small isolated clusters of correlated spins, and no collective relaxation 
associated with either dynamic or static transition is needed. 



V. THE GAUSSIAN POTTS GLASS 

We now turn our attention to the case where the couplings between neighboring 
spins are not given by a bimodal distribution, but by a Gaussian one (see Eq. (^)). 
As we will see, this change in the distribution has important implications for the 
static and dynamic properties of the system. All results discussed in this section 
were obtained for L — 10. 

In Fig. [L4] we show the T— dependence of the energy per spin. In agreement 
with our finding for the ±J model, this function does not show any sign for the 
presence of a singularity. At low temperatures we see that (with decreasing T) the 
slope of the curve decreases and that e(T) approaches smoothly its value in the 



ground state. This value is close to -5.49, the energy of the system if one assumes 
that all interctions are satisfied, i.e. 

= 3 (J )./,//•/,, - (p - 1) PWadJ^j . (23) 

Here P{Jij) is the distribution given by Eq. (0). The value of e^ pT is marked 
in the figure by the horizontal arrow. 



That T— dependence of the specific heat is shown in Fig. [jrj. Similar to the 
case of the ±J model, we find a Schottky-like peak at around T ~ 2. Since 
in the Gaussian model there is no gap in the excitation spectrum, by analogy 
with Ising spin glasses 0, we expect a low temperature behavior c(T) oc T and 
e(T)— e(0) oc T 2 rather than the behavior proportional to exp(— pJ/T) encountered 
in the ±J model. However, since for T < 1.2 the relaxation times are already 
rather large, no attempt has been made to verify these expectations numerically. 

Since the ferromagnetic bonds now occur with a broad spectrum of energies 
rather than a single energy pj, the argument put forward in the previous section 



for the existence of a plateau does no longer apply. In Fig. ^ we show the time- 
dependence of C(t) for all temperatures investigated and from this graph we see 
that in this T— range there is no sign for the presence of a plateau. Although for 
long times and low T the TTSP seems to hold reasonably well, we have not been 
able to determine the functional form of the time dependence of C(t). 

Of course it might be that the absence of the plateau is related to the fact that 
we are investigating the system at too high temperatures. In order to check for 
this possibility we have quenched 10 independent sample (i.e. different realisations 
of disorder) to T = 0.5 and have annealed the system for 10 7 MCS before we 
started a production run of 10 7 MCS. Although this time is certainly insufficient 
to equilibrate the systems at this temperature, it can be expected that for times 
significanly (O(10)) shorter than 10 7 MCS the correlation function is close to the 
one in true equilibrium. 

The correlation function as calculated after the quench is included in Fig. [ITJ 
as well. We see that even at this low temperature there is no well defined plateau. 
Instead the time dependence seems to be compatible with a logarithmic decay, at 
least in the time interval 10 < t < 10 5 . Hence we have evidence that for this model 
the mean field prediction for the existence of a dynamical transition does not seem 
to hold. 

If we define relaxation times similarly as done previously by requiring that C (t) 
reaches a particular value, 



C{t = n) = 0.1, C(t = t 2 ) = 0.05, C(t = r 3 ) = 0.02 (24) 

(horizontal dashed lines in Fig. [16]) we find again that the T; have an Arrhenius 
temperature dependence (see Fig. 4 in Ref. p8|), but the activation energy seems 
to increase smoothly the smaller the constant in Eq. (|24|) is chosen. (We find 
the activation energies 20.1, 16.4, and 12.8 for T3, T2 and t±.) Once more, this 



Id 



behavior is again expected from a system where one has small clusters (dimers, 
trimers, etc.) of ferro magnetically coupled Potts spins, and the breaking of this 
(continuous spectrum) of ferromagnetic bonds causes the extended tails of the 



time-correlation functions C(t) as seen in Fig. 16 



VI. CONCLUDING REMARKS 

In this paper we have used Monte Carlo simulations to investigate the static and 
dynamic properties of the 10-state Potts glass model with nearest neighbor random 
interactions on the simple cubic lattice. We have found that in the temperature 
range investigated there is no sign for the presence of the dynamical transition from 
ergodic to nonergodic behavior nor for the static discontinuous glass transition that 
one finds in the mean field version of this model. Instead the dynamic transition is 
replaced by a very gradual onset of slow dynamics, described by relaxation times 
that are compatible with simple Arrhenius laws. 

For the case of the ± J model, the spectrum of relaxation times is essentially 
discrete - namely times oc exp(kpJ/T) with k = 1, 2. These times are related to 
processes in which a spin is connected to its nearest neighbors via k ferromagnetic 
bonds that have to be broken for the spin to flip. In contrast to this, the Gaussian 
model has a spectrum of relaxation times that is continuous, which reflects the 
fact that also the values of the ferromagnetic bonds are distributed continuously. 

A further important difference between the ± J model and the Gaussian model 
is that the average spin auto-correlation function of the latter shows, after a fast 
drop at short time, a smooth decay to zero. In contrast to this, C(t) for the 
±J model shows several steps, the existence of which is related to the discrete 
spectrum of relaxation times. The auto-correlation function of individual spins 
is for most spins a simple exponential, with a decay time that is given by one of 
the relaxation times mentioned above. Thus, although there is a pronounced 
dynamical heterogeneity in the model, this heterogeneity is localized on a very 
small length scale (some spins have only antiferromagnetic bonds to their neighbors 
and can relax fully unhindered, others are in ferromagnetically coupled dimers, 
trimers, etc., but the concentration of large clusters is vanishingly small). Due 
to the strong localization of these excitations it is unlikely that a description of 
dynamical heterogeneity in terms of mesoscopically large regions ( "droplets" where 
volume and surface terms in their free energy compete, etc. [^, |37|) is applicable, 
at least for the present model. Instead we have shown in the present work that 
a description in terms of small clusters of ferromagnetically coupled spins in a 
background of spins to which they are "coupled" antiferromagnetically is almost 
quantitatively accurate. (Note that this antiferromagnetic coupling has almost 
no effect on the cluster, due to the large number of states for the spins that 
still are available when a bond between them is antiferromagnetic.) E.g. this 
description allows to explain the T— dependence of static quantities like energy, 
or the specific heat very well, and thus we conclude that in the present model 



the spins are correlated only very weakly The fact that we are able to give an 
almost quantitative description of the static as well as dynamic properties of the 
system in the T— range considered makes it also very unlikely that a dynamical 
or static transition occur at even lower temperatures, although strictly speaking 
such a possibility can of course not be ruled out. 

Hence we conclude that the idea that short range Potts glasses may provide 
a good model for the structural glass transition problem, exhibiting a rounded 
version of a dynamic transition (with a crossover from a power-law divergence 
of the relaxation time to an Arrhenius behavior with a very large barrier) and a 
(very closely avoided) static transition to a glass phase at a lower but nonzero 
temperature can not be maintained, at least not for the present model. Instead we 
think that a better model for the structural glass transition would be a Potts glass 
with interactions that have a somewhat longer range than in the models studied 
here. 

Very recently it has been suggested that such a close analogy between Potts 
glasses and structural glasses might exist only if the number of available states 
is very large, e.g. p = 10 3 [j37fl . However, we do not see any physical basis for 
such a claim: The choice AJ = 1, J = —1 in the model, necessary to avoid the 
tendency for ferromagnetic ordering, leads to a fraction x ~ 0.146 of ferromagnetic 
bonds, and for larger values of p an even smaller x is needed [0. However, 



with such a small fraction of x the properties of the model will be qualitatively 
the same for p = 10 and p = 10 3 . Frustration effects also cannot change this 
conclusion - the fraction of frustrated plaquettes with three ferromagnetic bonds 
and one antiferromagnetic bond simply is too small (note that plaquettes with 
one ferromagnetic bond but three antiferromagnetic bonds are frustrated in an 
Ising spin glass, where p = 2, but not for the large p values considered here). 
While the idea |J7| that one should discuss the applicability of the mean field 
Potts glass transition scenario to real systems in terms of a Ginzburg criterion 



is certainly interesting, it seems that it is not useful for the extremely short 
range case of a nearest neighbor Potts glass. This conclusion is of course not 
that surprising since our experience with Ising ferromagnets has shown that the 
Ginzburg criterion tells little about the nearest neighbor model either, although it 
is useful for understanding properties of Ising models with a large but finite range 
of the forces ||53|| . A study of the 10— state Potts glass with random interactions 
with such medium range could be illuminating, but clearly is computationally 
very demanding and hence was not attempted here. In any case, our study shows 
that it is rather difficult to draw analogies between spin-glass-type models and 
models commonly used to study the problem of the structural glass transition 
(hard spheres, Lennard- Jones, etc.) even if the description of such models within 
mean-field looks very similar. 
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APPENDIX 



In this appendix we briefly describe the analytical calculations we used to obtain 
the (approximate) temperature dependence of the energy e(T) and of the specific 
heat c(T). 

To calculate e(T) we consider a Potts-spin that we will denote by "0". This 
spin is connected to 2d neighbors (d is the dimension of the lattice, i.e. in our case 
d — 3). Assume that out of the 2d bonds that connect the central spin to its 2d 
neighbors, exactly k are ferromagnetic and that hence 2d—k are anti-ferromagnetic. 
In the following we will calculate the partition function for such a "fc— cluster" of 
2d + 1 spins and subsequently its energy. We then make the approximation that 
the Potts Hamiltonian given by Eq. (|I]) is a sum of independent clusters (using the 
distribution P{Jij) to determine the frequency with which a k— cluster occurs). 

Using the notation A := pj, we can write the Hamiltonian of a A;— cluster as 

k 2d 

Hk = -A ^2 <W< + A Y. 6 ^ + 2 J ( k ~ d ) ■ ( Al ) 

i=l i=k+l 

The partition function can now easily be obtained as 

Z k =p[p-l + exp(/3A)] fc [p - 1 + exp(-/3A)] 2d ' k exp[(32J(k - d)} , (A.2) 

with (3 = 1/ksT. The energy can now be calculated by using the relation efc(T) = 
—Z^dZ/dp and one obtains: 



et(T) 



A f (M-t)exp(-M) kexp(PA) \ 

\ p - 1 + exp(-^A) p - 1 + exp(/3A) J 



The specific heat ck(T) can now be calculated using c^(T) = dexiT) / dT which 
nves 



c k (T) = A 2 (3 2 (p-l)exp(-(3A) 



(2d - k) 



[p-l + exp(-/3A)] 2 {(p - 1) exp(-/3A) + l]) 2 

(A.l) 

If in the real system we pick an arbitrary spin and consider its 2d nearest 
neighbors there is a probability W(k) that the first spin has exactly k ferromagnet 
bonds and W(k) is given by: 



(A.5) 



1U 



where x is the concentration of ferromagnetic bonds (see Eqs. (0) and (^)). Thus 
within the approximation considered here, the value of an observable Y of the 
system (i.e. energy per spin or specific heat) is given by 



2d 



Y = Y,YuW{k) . (A.6) 



=o 



Using this relation we have calculated the temperature dependence of the energy 
per spin as well as the specific heat that are shown in Figs. ^ and respectively. 
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FIG. 1: Schematic sketch of the thermodynamic and dynamic transitions in a mean- 
field p-state Potts glass model with p > 4 states. The spin glass order parameter qEA 
is non-zero only for T < Tq and jumps to zero discontinuously at T = Tq. The spin 
glass susceptibility xsg follows a Curie- Weiss type relation with an apparent divergence 
at a "spinodal temperature" at T s < Tq, and hence is still finite at To. The relaxation 
time r diverges according to a power-law r oc (T — To) _A at a temperature To > To. 
This divergence is due to the occurrence of a long-lived plateau of height qEA in the 
time-dependent spin auto-correlation function C(t). The entropy s(T) per spin has no 
singularity at To, but shows a kink singularity at To (thus there occurs no latent heat, 
which would mean an entropy jump). The extrapolation of the high temperature branch 
of the entropy would vanish at a temperature Tk < Tq. 
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FIG. 2: Temperature dependence of the ferromagnetic susceptibility xfm = ^V[("i 2 )]av 
in the p = 10 Potts glass for (a) Jo = and (b) Jo = — 1. The three curves correspond 
to different system sizes. 
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FIG. 3: Temperature dependence of the energy per spin e(T) for three lattice sizes 
as indicated in the figure. Note that for T > 1.5 the average [■ • • ] av was realized by 
averaging over 100 independent realizations for L = 6 and L = 10, and by averaging 
over 50 independent realizations for L = 16. For T = 1.3 and T = 1.4, only L = 10 was 
studied, and only 10 realizations were used, since 10 times longer runs than usual were 
necessary. The solid curve is the theoretical prediction using an approach of independent 
clusters (see main text for details). 




FIG. 4: Specific heat per spin c(T) of the bimodal Potts glass plotted vs. temperature for 
three lattice sizes, as indicated in the figure. The bold curve is the result of a theoretical 
calculation using an approach of independent clusters (see main text for details). The 
inset shows a lin-log plot of c(T) vs. 1/T, in order to show the T— dependence at low 
temperatures. 
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FIG. 5: Scaled first moment N 1 / 2 ^) of the order parameter plotted as a function of 
temperature. Different curves correspond to the different choices of L, as indicated. The 
inset shows the T— dependence of (q). 
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FIG. 6: Spin glass susceptibility xsg as defined in Eq. (|l3|) as a function of temperature. 
The inset shows the same data in a log- log plot, to demonstrate that at high T we have 
Xsg oc 1/T. 
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FIG. 7: Plot of the reduced fourth-order cumulant gn{N,T) versus temperature, for 
L = 10 (open circles) and L = 16 (filled squares). The data for L = 6 were rather noisy 
and therefore were not included. 
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FIG. 8: Plot of the Guerra parameter G(N, T) versus temperature for L = 10 (open 
circles) and L = 16 (filled squares). 




t [MCS] 

FIG. 9: Time dependence of the spin auto-correlation function of the ± J Potts glass for 
L = 10 and various temperatures (left to right: T = 2.4, 2.2, 2.0, 1.8, 1.7, 1.6, 1.5, 1.4, and 
1.3). The three horizontal dashed lines are used to define the relaxation times Tj (see 
main text for details). In the inset we compare C(t) for L = 10 and L = 16 at T = 1.6. 




FIG. 10: a) Spin auto-correlation function C(t) for L = 10 plotted vs. t/r\. b) Spin 
auto-correlation function C(t) for L = 10 plotted vs. t/r2- The data for C(t) of this 
figure are identical to the ones shown in Fig. |||. 
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FIG. 11: Arrhenius plot of the relaxation times t\, T2, and t% defined in Eq. (|l9|), for the 
± J Potts glass. The straight lines are fits with an Arrhenius law and from their slopes 
the activation energies E% , see Eq. (20), are estimated. 



Z,<J 




t [MCS] 

FIG. 12: Time dependence of the relaxation function Cj(i) = (p— l)~ 1 (S'i(t / ) • Si(t' + t)) 
of individual spins of the ± J Potts glass for L = 10 at three different temperatures: 
(a) T = 2.4, (b) T = 2.0, and (c) T = 1.6. 
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FIG. 13: Normalized distribution of barrier heights bi and pre-factors (inset) obtained 
from fitting the Arrhenius relation, Eq. (pi]), to the relaxation times t[ of individual spins, 
using either the two temperatures T = 1.6,2.0 (full curves) or T = 2.0,2,4 (dashed 
curves) . 
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FIG. 14: Energy per spin e(T) of the 10-state Gaussian Potts glass plotted versus 
temperature for L = 10. The arrow indicates the estimate of the ground state energy 
under the assumption that there are no frustration effects (see Eq. (|2l 
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FIG. 15: Specific heat c(T) per spin of the 10— state Gaussian Potts glass plotted versus 
temperature for L = 10. Points are from energy fluctuations, while the broken line is a 



numerical derivative of the curve e(T) shown in Fig. 14 
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FIG. 16: Spin auto-correlation function C(t) of the Gaussian Potts glass plotted versus 
time for L = 10. The temperatures are (from left to right): T = 2.8, 2.6, 2.4, 2.2, 2.0, 
1.8, 1.7, 1.6, 1.5, 1.4, 1.3 and 0.5. The horizontal dashed lines are used to define the 
relaxation times Ti given in Eq. (p4]). Note that the data for T = 0.5 is not in equilibrium 
(see main text for details). 



